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We explore the photonic (bright) side of dynamical Coulomb blockade (DCB) by measuring the 
radiation emitted by a dc voltage-biased Josephson junction embedded in a microwave resonator. 
In this regime Cooper pair tunneling is inelastic and associated to the transfer of an energy 2eV 
into the resonator modes. We have measured simultaneously the Cooper pair current and the 
photon emission rate at the resonance frequency of the resonator. Our results show two regimes, in 
which each tunneling Cooper pair emits either one or two photons into the resonator. The spectral 
properties of the emitted radiation are accounted for by an extension to DCB theory. 

PACS numbers: 74.50+r, 73.23Hk, 85.25Cp 



Dynamical Coulomb blockade (DCB) of tunneling is 
a quantum phenomenon in which tunneling of charge 
through a small tunnel junction is modified by its elec- 
tromagnetic environment [IHl]. This environment is de- 
scribed as an impedance in series with the tunnel element 
(see Fig. [T^) . The sudden charge transfer associated with 
tunneling can generate photons in the electromagnetic 
modes of the environment. In a normal metal tunnel 
junction, biased at voltage V, the energy eV of a tunnel- 
ing electron can be dissipated both into quasiparticle ex- 
citations in the electrodes and into photons. At low tem- 
perature energy conservation forbids tunneling processes 
emitting photons with total energy higher than eV . This 
suppression reduces the conductance at low bias voltage 
[TJI21I1]- In a Josephson junction, DCB effects are more 
prominent since at bias voltages smaller than the gap 
voltage 2A/e quasiparticle excitations cannot take away 
energy. Therefore, as shown in Fig. [T^, the entire energy 
2eV of tunneling Cooper pairs has to be transformed into 
photons in the impedance for a dc current to flow through 
the junction [3 E] . Experiments have confirmed the pre- 
dictions of DCB theory for the tunneling current, both 
in the normal [IHZ] and superconducting case [3 |9] but 
the associated emission of photons into the environment 
has never been investigated. The aim of this work is pre- 
cisely to fill this gap by exploring the photonic side of 
DCB. We do so by embedding a Josephson junction into 
a well controlled electromagnetic environment provided 
by a microwave resonator. The resonator in turn leaks 
photons into an amplifier, allowing to measure the rate 
and spectrum of photons emitted by the junction. 

The experimental setup is represented in Fig. [T|d. A 
small SQUID acts as a tunable Josephson junction with 
Josephson energy E'j — E'jqI cos(e<I>//i)| adjustable via 
the magnetic flux $ threading its loop. The microwave 
resonator is made of two quarter-wave transformers and 
its fundamental mode has frequency ~ 6.0 GHz and 
quality factor Qo — 9.4. Higher modes of the resonator 
appear at i/„ ~ (2n -I- l)z^o {n = 1,2,.. .) with the same 
lineshape up to small deviations caused by the junction 
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FIG. 1: Dynamical Coulomb blockade of Josephson tuimcl- 
ing. a-Generic circuit: a dc voltage-biased Josephson element 
in series with an electromagnetic environment with impedance 
Z{p), which can be described as a collection of harmonic os- 
cillators. Below the gap voltage, the energy 2eV of a tun- 
neling Cooper pair is transferred to the electromagnetic en- 
vironment in form of one or several photons in its modes. 
b-Experimental setup: The sample consists of a SQUID, 
working as a tunable Josephson junction, in series with a mi- 
crowave resonator, forming the electromagnetic environment 
and consisting of two quarter-wave transformers. A bias tee 
separates the low frequency voltage bias V and current mea- 
surement / from the emitted microwave power P^jw Three 
circulators (only one represented here) ensure thermalization 
of the microwave environment. The emitted power is ampli- 
fied at 4.2 K and then bandpass filtered and detected at room 
temperature. 



capacitance (estimated to be 4fF). The expected total 
impedance Z{v) seen by the Josephson element is plotted 
in the top panels of Figs. [2] and [3) and reaches ~ 1.5 kfi 
for all modes. 
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The sample is cooled to T ~ 15 mK in a dilution re- 
frigerator and connected to a bias tee separating the 
high-frequency and low-frequency components of the cur- 
rent. The low- frequency port is connected to a voltage 
bias through a Ikfi resistor used to measure the tun- 
neling current /. The microwave port is connected to a 
4 to 8 GHz cryogenic amplifier with noise temperature 
Tn — 3.5 K through three circulators protecting the sam- 
ple from the amplifier noise and ensuring thermalization 
of the environment. After further amplification and even- 
tually bandpass filtering at room temperature, the power 
is detected using a calibrated square-law detector. Its 
output voltage is proportional to the microwave power, 
which contains the weak sample contribution P^^ on top 
of the large noise floor of the cryogenic amplifier. In or- 
der to measure P^w, we remove the background with a 
low frequency lock-in technique by applying a to V 
square-wave bias voltage modulation to the sample. 

We first characterize the on-chip microwave resonator 
and determine the gain of the microwave chain by mea- 
suring the power emitted by the electronic shot noise of 
the junction Sn ~ el at bias voltage V ~ 2mV, well 
above the gap voltage 2A/e ~ 0.4 mV. Under these 
conditions, the spectral density of the emitted power is 
2eVRe Z{iy)Rj^/\Rj^ + Z{iy)\^ ~ 2eV Re Z (v) / with 
i?N = 17.9 kfi ^ \Z(v)\ the tunnel resistance of the junc- 
tion in the normal state. This spectral density is ex- 
tracted using a heterodyne measurement equivalent to 
a 100 MHz wide bandpass at tunable frequency. More 
details about this procedure can be found in |10| . The 
extracted Re Z(y) is shown in the upper panel of Fig. [2] 
and the right panel of Fig. [4] It shows the expected peak 
at Vq with an additional modulation due to parasitic re- 
fiections in the microwave set-up. 

Once the system calibrated, we measure the Cooper 
pair transfer rate across the SQUID F'-'p — I /2e and the 
photon emission rate into the fundamental mode of the 
resonator Fq*^ as a function of bias voltage. Rate Fq*^ 
is extracted from the microwave power P^wi emitted in 
a 2 GHz wide band centered at va = 6 GHz, via Fq*^ = 
Pfivr/hvo- In order to simplify interpretation of our data, 
we adjust iSj for each experiment to ensure good signal 
to noise ratio while keeping the maximum Fq'' smaller 
than the decay rate of the resonator 2'kvq/Qq ~ 4 GHz, 
so that the electromagnetic environment of the junction 
stays close to thermal equilibrium. 

The measured rates F'-'p and Fq'' are shown in Figs. [2] 
and[3j In Fig.[2]both rates display a peak at voltage Vq ~ 
hvQ/2e = 12 /iV, where the energy of a tunneling Cooper 
pair corresponds to the energy of one photon at frequencv 
vq. On resonance F^p and Fq^ agree within 5 %. In Fig.p^ 
additional peaks appear in F'-'p at voltages Vn = /it'„/2e 
(n = 1,2,3), which do not have counterparts in Fq^. In 
both figures F^p and Fq'' show much smaller peaks at 
V^ + Vr, (n = 0,l,2,3). 
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FIG. 2: Fundamental mode uo of the resonator and corre- 
sponding Cooper pair and photon rates. Top panel: real part 
of the impedance seen by the junction, calculated from the 
resonator geometry (black line) and reconstructed (magenta) 
from a quasiparticle shot noise measurement [10]. Bottom 
panel: measured Cooper pair rate F^^p (red) and photon rate 
Fq*^ (blue) integrated from 5 to 7 GHz. Both rates show a 
peak around V = hvo /2e ~ 12 ^V. A second small peak in the 
microwave power at V — 24 fiV (vertical zoom in inset) corre- 
sponds to two-photon processes. Solid lines are P{E) theory 
fits using the reconstructed impedance from 5 to 7 GHz and 
the calculated one outside this range, an effective temperature 
of 60 mK, and a single adjustable parameter _Bj = 5.1 /ieV. 
Cyan and yellow hues correspond, respectively, to Eq. ([T]) and 
Eq. ||5| integrated from 5 to 7 GHz. 



DCB theory accounts for these results. At low tem- 
perature kgT <^ 2eV , when backward tunneling is neg- 
ligible, the Cooper pair rate is given by OH]: 



= ^E]P{2eV), 



(1) 



where P{E) is the probability density for a tunneling 
charge 2e to emit an energy E in form of photons into the 
impedance. The function P{E) is the Fourier transform 
of exp J{t) with |4j 



J{t) ^2 I dv 



r[v) exp(— i27ri/i) — 1 
1 — exp(— /ii^/fceT) 



where 



r{v) 



4e^ 



ReZ{u). 



(2) 



(3) 



At T = 0, the expansion of the function P{E) in powers 
of r{iy) yields: 
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FIG. 3: Higher orders modes of the resonator and corre- 
sponding Cooper pair and photon rates. The plotted quan- 
tities are the same as in Fig. [2j but for a larger voltage 
and frequency range and different Josephson coupling. The 
Cooper pair rate shows peaks at voltages V — hvn/2e, as- 
sociated to the emission of photons in the resonator mode 
at frequency i^n- The power measurement only probes the 
photons emitted in the fundamental mode of the resonator. 
It shows a large peak a.t V = hi'o/2e, and smaller ones at 
V — h^vo + Vn) I'ie, associated to the simultaneous emission of 
one photon at and one photon at u„ per tunneling Cooper 
pair. Here theoretical curves are obtained from the calcu- 
lated impedance in the whole frequency range and the fitting 
parameter Ej = 10.3 ^eV. 



where = (I - Jdur{v)/v) and ^ 2eV/h. We 
neglect here the low-frequency environment and suppose 
^dvr{v)/v <C 1. The first term in the parenthesis of 
Eq. Q corresponds to the emission of a single photon 
at frequency vj. It accounts for the peaks in F^p ob- 
served at bias voltages Vn, and for the peak in Fg*^ at Vq. 
Since we only detect photons emitted in the fundamen- 
tal mode, Fq*^ does not show any peak at 14 for n > 0. 
The second term in the parenthesis corresponds to the 
simultaneous emission of two photons at frequencies v 
and i^j — v. It accounts for the peaks in Fq'' observed 
at bias voltage + Vb, where the energy of a tunneling 
Cooper pair can be split into one photon at frequency 
and another at frequency Vn- These second order pro- 
cesses have much lower rates than the first order peaks 
in F'-'p because of the low resonator impedance. The 
same processes also appear in F^p as small peaks be- 
tween the main resonances. Note however that first order 
processes contribute about equally to the peaks in F^p 
due to small rises in the impedance at vq -(- Vn. Eq. |4] 
also accounts for the increase in F^p at low voltage in 
Fig. [2] and the reduction of the height of the resonances 
at Vn with increasing n, which both reflect the 1/v scal- 
ing of the first order term. Note that as far as single 
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FIG. 4: Main panel: spectral density of the photon emission 
rate 7 for the first and second order peaks {i/q and 2i/o) as a 
function of bias voltage and frequency. Data were taken for 
difi'erent Josephson energies and were normalized by the max- 
imum 7 of each data set. Dotted lines correspond to 2eV = hi/ 
and 2eV = h{i' + vq). White dots correspond to power spec- 
tra at fixed voltages (12.5 and 24.5 /xeV) and solid lines to 
theoretical predictions Eq. Q integrated over 100 MHz, the 
width of the filter used in this experiment with an effective 
temperature of 60 mK and fitting parameters Ej = 5.0 fieV at 
V = 12.5 and Ej = 14.3 fieV atV = 24.5 fiY. Right panel: 
reconstructed dissipative part of the resonator impedance, as 
in Fig. (2] 

photon processes are concerned, the ac Josephson effect 
predicts the same voltage to frequency conversion and 
Cooper pair and photon rates. However, DCB assumes 
that the transferred charge is almost a classical variable 
and its conjugate variable, the phase across the junction, 
has large quantum fluctuations, whereas the opposite as- 
sumption is made in the ac Josephson effect. 

In Fig. [4] we explore the spectral properties of the 
emitted radiation, using a heterodyne technique with 
100 MHz resolution. The first order emission occurs along 
a narrow line at 2eV = hv whose amplitude is modu- 
lated by r{v), in agreement with Eq. This equation, 
however, makes no prediction about the spectral width of 
the emitted radiation which is found to be approximately 
150 MHz, significantly narrower than the 700 MHz wide 
resonance in r(i^). The second order peak follows a 
2eV = h(iy + vq) line and has a much broader spectral 
width, comparable to the width of the resonance in r(^'). 

In order to account for these observations we calculate 
the spectral density of photon emission in the DCB 
regime, which is related to the emission current noise 
spectral density Su via 7 = 2'ReZ{y)Sii{v)/hv. From 
the relation between Su and P{E) [31 [T^ one obtains 

^^^I^.^EjP{2eV-h.). (5) 
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This expression, valid at low temperature fc^T ^ 
hiy, 2eV, can be understood within a Caldeira-Leggett 
decomposition [11] of the impedance Z{i^) into modes 
with infinitesimal width |10| : The term 2r{v)/i' describes 
the probability density for each tunneling Cooper pair to 
emit a photon at frequency v into the environment. The 
remainder of Eq. ([s]) gives the Cooper pair rate while 
emitting the remaining energy 2eV — hv in the form of 
one or several other photons, as in Eq. ([!]). 

The first order peak occurs at 2eV ~ hv^ so that P{E) 
has to be evaluated around i5 ~ where the perturba- 
tive expression Q is not valid. Therefore we calculate 
P{E) at all orders in r at finite temperature, following 
the method described in Ref. [H]. P{E ~ 0) is approx- 
imately a Lorentzian jyii^^^ of width 2D ~ 47rfcTr(0) 
. This expression gives the spectral width of the first 
order peak. Its height along the line 2eV = hv follows 
the variations of r(v) due to the first term in Eq. ([5]). 
In order to find a quantitative agreement with the data 
we use an effective temperature T = 60 mK, in good 
agreement with the 60 mK ± 10 mK electronic tempera- 
ture deduced from shot noise measurements performed 
at vq on the same junction, driven to the normal state 
by applying a 0.1 T magnetic field [13]. We attribute the 
difference between the effective temperature of the envi- 
ronment and the fridge temperature of 15 mK mainly to 
low frequency noise radiated by the amplifier used for the 
current measurement. 

For the second order peak P{E) has to be evaluated 
around E ~ hvQ, where it reproduces the variations of 
r{E/h). The radiation of the second order peak is there- 
fore emitted in a band around 2eV = h{vo + ^) with a 
spectral width given by the width of the resonance in 
r(j/). As for the first-order peak, the amplitude along 
this line is then modulated by r(y) due to the first term 
in Eq. ^. 

With the spectral distribution 7 at hand we now come 
back to the data presented in Figs. [2] and [3] for a more 
quantitative description. The solid yellow lines in both 
figures show expression ([5]), integrated from 5 to 7 GHz, 
and the solid cyan lines expression 

US- The theory 
accurately matches the measurement with a single fitting 
parameter, the Josephson energy and accounts for 
the relative weight of the first and second order peaks in 
Fig. [2j Discrepancies only show up at the larger bias volt- 
ages in Fig. [3] corresponding to frequencies above 18 GHz 
where we do not accurately control the electromagnetic 
environment. For the sake of consistency with Fig. |4] we 
have used the non perturbative P{E), although the zero 
temperature perturbative expression Q would already 
have given a satisfactory agreement with the data. 

In conclusion, our results clearly validate the under- 
standing of DCB based on photon emission into the elec- 
tromagnetic environment. We have measured the spec- 
tral density of the emitted radiation, which we quantita- 
tively account for by an extension to DCB theory. A nat- 



ural continuation of this work is the investigation of tun- 
neling current and emitted radiation at higher Josephson 
energy, where the emission rate can be strong enough to 
overcome cavity loss and cause stimulated emission. A 
larger Josephson energy will also make simultaneous tun- 
neling of multiple Cooper pairs relevant. These are the 
two ingredients for the transition from the DCB regime 
to the regime of the ac Josephson effect. Furthermore, 
our results demonstrate that under appropriate biasing, 
the tunneling of a Cooper pair is associated with the si- 
multaneous emission of multiple photons. This opens the 
way to new schemes for the emerging field of microwave 
quantum optics using inelastic tunneling of Cooper pairs 
for the generation of non classical photon states. 
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IdF, technical help from P. Senat, P. Jacques, P. Orfila 
and discussions within the Quantronics group. 
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SUPPLEMENTARY MATERIAL 
Details on setup and calibration 

We describe the calibration of the low frequency cir- 
cuitry for voltage bias and current measurement and the 
microwave components used to define the environment of 
the junction and to measure the emitted radiation. 



Low frequency circuit 

In addition to the components depicted in Fig. 1 in the 
paper, the low frequency circuit includes multipole RC 
low pass filters in the bias line (cutoff at 8 kHz) and the 
current measurement lines (cutoff at ~ 100 kHz), several 
capacitors to ground and a copper powder filter between 
readout resistor and bias T (see Fig.js]). 

Calibration of the division factor of the voltage di- 
vider used for biasing the sample is performed in situ by 
reading out the applied bias voltage through one of the 
lines used for current measurements (see Fig. 1). We use 
a lock-in technique for this measurement and minimize 
the critical current of the junction so that the current 
through the sample is negligible. 

The 1 kil NiCr resistor used for current measurement 
cannot be calibrated in situ due to large lead resistances, 
but we have verified that a resistor from the same batch 
does not deviate from its nominal value by more than 
1 % at the operating temperature of our refrigerator. 

In the presented measurements we have corrected the 
applied bias voltage for the voltage drop over the 1 kil 
resistor. 

The good agreement between theory and data in Fig. 4 
supports the accuracy of this DC calibration. 



Microwave circuit 

The microwave chain (see Fig. [5| consists of the junc- 
tion, the resonator, a bias T, a 6dB attenuator, a 
triplexer with crossover frequencies 4 and 8 GHz, two 4 
to 8 GHz isolators at base temperature, an isolator at 
700 mK, a cryogenic amplifier at 4K with 3.5 K noise 
temperature, a 3 dB attenuator at 4 K followed by a room 
temperature amplification and measurement chain. The 
isolators protect the sample from the amplifier noise in 
their working range from 4 to 8 GHz. The triplexer, work- 
ing from DC to 18 GHz, connects the sample to the mi- 
crowave chain only in the 4 to 8 GHz range. Outside 
this range the sample is connected to cold 50 ft loads so 
that the sample "sees" a cold 50 load throughout the 
18 GHz frequency range of the triplexer. Refiections due 
to the insertion loss of the triplexer are suppressed by the 
6dB attenuator. 
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FIG. 5: More detailled view of the experimental setup. Only 
circuit components inside the refrigerator are shown in full 
detail. 



Calibration of the resonator impedance Z{u) and the 
gain of the microwave chain from the sample to the top of 
the refrigerator is performed in situ because the super- 
conducting resonator must be cooled below its critical 
temperature and the amplifier gain and cable attenua- 
tion depend on temperature. We perform both calibra- 
tions using quasiparticle shot noise as white current noise 
reference. We bias the junction at 2mV, far above the 
gap voltage 0.4 mV of our junction, where the current 
noise is in good approximation Su = el at frequencies 
^ eV/h ~ 0.5 THz. In order to separate this noise 
from the noise fioor of the cryogenic amplifier, we then 
apply small variations of the bias voltage and measure 
the corresponding changes in the measured microwave 
power with a lock-in amplifier. 

The conversion of Su into emitted microwave power 
depends on the environment impedance Z{i') seen by the 
junction and its tunneling resistance i?N- First, only a 
fraction i?^/|i?N + ^(i^)P of the current noise is absorbed 
by the environment, but this factor is close to 1 (> 0.85) 
as i?N = 17.9 krj is much larger than the environment 
impedance. The current noise in the environment has 
then to be multiplied by 2ReZ(i/) to obtain the mi- 
crowave power emitted by quasiparticle shot noise: 



Sp{iy) = 2eV 



ReZ(^y)i^^ 



12 ■ 



(6) 



The emitted power depends on Zlv), therefore Z(iy) and 
the gain of the microwave chain cannot be calibrated in- 
dependently. We solve this problem by making the fol- 
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lowing assumptions: 

1. The gain of the microwave chain is frequency inde- 
pendent in the 5 to 7 GHz range. 

2. The integral over the peak in Re Z^i/) from 5 to 
7 GHz corresponds to calculations based on the res- 
onator geometry. 

Assumption 1 is justified by a previous calibration of the 
gain of the amplifiers, found to be constant in the 5 to 
7 GHz range within 0.5 dB, and by observing that the 
noise floor of the amplifier in this experiment in the same 
5 to 7 GHz range deviates by less than 0.5 dB from its 
mean value, meaning that the gain and attenuation af- 
ter the input of the cold amplifier are almost frequency 
independent. However, frequency dependent attenuation 
between the sample and the cryogenic amplifier is not 
detected this way and can be source of slight inconsis- 
tencies. 

Assumption 2 relates to the fact that the integral over 
ReZ(i/), proportional to the characteristic impedance 
of the resonance, mainly depends on the transmission 
line impedance of the quarterwave segment closest to 
the junction, less on the second segment, and only very 
weakly on the load impedance beyond the resonator. For 
example, the characteristic impedance of the resonance 
changes by less than 20 % when the the load impedance is 
changed from 50 57 to any value between 20 f2 and infinity. 
Therefore we can get an accurate estimate of this inte- 
gral from the impedance of the resonator segments, even 
in the presence of imperfections in the line impedance 
seen from the resonator. We calculate the character- 
istic impedance of our quarterwave segments from the 
designed resonator geometry, checked with an electron 
microscope and the well known stackup of our sample. 
We obtain characteristic impedances 13917 and 25.5 51, 
resulting in a characteristic impedance of the resonance 
mode of 150 il. The fact that P{E) theory using this res- 
onator impedance accurately fits the relative heights of 
first and second order peaks in Fig. 2 indicates that this 
estimation of the resonator impedance is correct. 

With these two assumptions we can use Eq. ^ to sep- 
arately calibrate the power gain and Ziy). To calibrate 
the power gain in Figs. 2 and 3 we compare the measured 
shot noise signal to Eq. ([6]) integrated over the bandwidth 
of our measurement ranging from 5 to 7 GHz. 

In order to determine the spectral dependence of Z[v) 
we use a narrow band pass filter (as for the data in 
Fig. 4) to measure a signal proportional to Sp. To ex- 
tract Ik^Ziy) from the emitted power using Eq. ([6]), one 
would need to know lmZ(v), which is not measured in- 
dependently. To circumvent this difficulty, we use Eq. ^ 
iteratively to calculate Z{i'): we start with the calculated 
impedance Re Z^iy) in the denominator and in each step 
calculate a new YleZniv), normalize it using assumption 
2, and then calculate Im Zn{v) using the Kramers-Kronig 



relations with ReZ„(z/) in the 5 to 7 GHz range where 
Sp was measured and ReZo(z^) outside this range. This 
procedure converges quickly and yields a precise determi- 
nation of Re Ziy) because Im Z{v) ^ and therefore 
only has a small influence on the extracted ReZ(j/). 



Noise emitted by a Josephson junction in the DCB 
regime 

We give here two derivations of expression (5) in the 
paper. The first one follows the derivation of Coulomb 
blockade in 4J and treats photon emission as fluctuations 
in the Cooper pair tunneling rate. The second derivation 
explicitly evaluates photon emission and absorption rates 
in a mode of the electromagnetic environment. 



Derivation 1 

The Josephson element is described by the Hamilto- 
nian ifj = — iJjcosJ and the current operator through 
the element is 



The current-current correlator function is then 



Suit) - (mm 



To arrive at the last line we have noted that at non-zero 
dc bias voltage V terms e^"'*^"-' with the same sign 

in the exponential average to zero. We decompose the 
phase difference across the tunnel element 5{t) = (5o(i) -I- 
5{t) into the deterministic part 5o{t) = 2eVt/h and a 
fluctuating random phase 5{t) caused by the fluctuations 
in the electromagnetic environment. We use the gaussian 
character of the noise of the linear environment to rewrite 
14]: 

g±M(t)gT*5(0)\^ ^ g±j2eVt/n ^g±j5(t)gT«(0)^ 
^ ^±i2eVt/h^J{t) 

where J{t) = (^{6{t) - ~6{0)^ ~5{Q)^ is the (supercon- 
ducting) phase-phase correlation function across the 
impedance, so that ([t]) becomes 

Suit) = -^2cos— e^W. 

The (non-symmetrized) current noise density Sii{i') is 
directly the Fourier transform of the current -current cor- 
relator Sii{t) 



Sii{u) 



Sn{t)e- 
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(Wiener Khinchin theorem) which yields immediately 



Siiiu, V) 



{P{2eV - hv) + P{~2eV - hv)) 



with P{E) = 2^ ^J{t)+iEt/h^^ ^Yie Fourier trans- 
form of exp J{t). Assuming that the mode of the environ- 
ment at frequency v is in the ground state, i.e. kT <Si hv, 
the current noise gives rise to a photon emission rate 
density 

2ReZH 
7 = ; '^11 



hv 



V 2h 



E] [P{2eV - hv) + P{-2eV - hv)) 



In expression (5) in the paper we have neglected the 
second term corresponding to photon emission during 
backward tunneling which is negligible for kT <C 2eV + 
hv as in our experiment. 



where p/^„ is the probability for mode I to be in state |n) 
before the tunneling event. We assume that the tunnel- 
ing rate is low so that the environment can reach thermal 
equilibrium before each tunneling event and the pi^n de- 
scribe thermal equilibrium states. This equation can be 
seen as the definition of P{E), the probability density 
for the environment to absorb energy E when a Cooper 
pair tunnels. Calculating P{E) by evaluating this sum is 
possible but cumbersome, so usually the approach from 
Ref. [Ij is taken. 

In order to calculate the rate for emitting photons at 
frequency v we treat mode k at Vk = v separately and 
keep track of the photon number difference ruk- We ab- 
sorb all other modes into a function P'{E) describing a 
modified environment where mode k has been removed. 



E 



Pk,n 
n— max{0. — m} 



n+m 



2h 



E]P'{2eV-mhvk). 



Derivation 2 
The Hamiltonian of the total system is 

H = Hera + H} 



with 



Hem — ^ 



{Ml 

2L, 



Ql 
2Ci 



Hj = -Ej cos S 
2eVt 



Here Hem is the Hamiltonian describing the environment 
of impedance Z(uj). Ci and L; describe the mode at 
frequency vi = 1/2tt\/CiLi. A second constraint on C; 
and Li imposed by the requirement that the sum of the 
resonator impedances approaches Z(v), giving 



Z, = 



Ci n V 



where Av is the spacing between adjacent modes, sup- 
posed to be constant. The Cooper pair tunneling rate 
can now be calculated using Fermi's golden rule: 



TT 

2h 



n; > I 
rii -{-mi > 









(ni + mi 







S^^mihvi - 2eVj 



2h 



EjP{2eV), 



In this expression the phase operator (pk across the res- 
onator k can be expressed in terms of the rising and low- 
ering operators a\ and ak'- 



with 



h 



4e2 ^ Av 

2—ReZ{v) 

h V 



2r{v) 



Av 



The spectral density of photon emission/absorption is 



where e*'^'= 



1 



calculated in the limit Av 
«V^(4 + «fe) and P'{E) P{E). In this limit r;;^ 
is nonzero only for m — 0,±1. The rate without pho- 
ton emission or absorption F,^ simply tends to the bare 



Cooper pair tunneling rate, i.e. Fg 



For ra = ±1 



we calculate the density of photon absorption/emission 



= 2 



lim — ^ 
Av^O Av 

CX) 



(8) 



r{v) 



n=l 



.-l{^)i^EiP{2eVThv). 



The backward tunneling processes F'*" and ^'^^ have the 
same expressions as the corresponding forward processes 
but with V replaced by — y. 

At low temperature kT <^ hv the mode at v is in the 
thermal ground state, i.e. Pn{v) ~ 5n,o, and P{—2eV — 
hv) ~ 0, so that 7!^ = 7^7 = 0. We then find expression 
(5) in the paper: 



2r{v) TT 
V 2h 



ElP{2eV - hv) 



